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, The Casimir energy is the first-order-in-/l correction to the energy of a 

time-independent field configuration in a quantum field theory. We study 
the Casimir energy in a toy model, where the classical field is replaced by 
a separable potential. In this model the exact answer is trivial to compute, 
■ making it a good place to examine subtleties of the problem. We construct 

two traditional representations of the Casimir energy, one from the Greens 
function, the other from the phase shifts, and apply them to this case. We 
show that the two representations are correct and equivalent in this model. 
We study the convergence of the Born approximation to the Casimir energy 
t>> and relate our findings to computational issues that arise in more realistic 

^ models. 

On 

I. INTRODUCTION 

Qh' 

Casimir energies are important in quantum field theories, where they give the first- 



order (in h) corrections to the energy of time-independent classical field configurations. 
Computations of the Casimir energy involve formal manipulations of divergent expressions 
that eventually are regulated and renormalized. [[| Since they go beyond perturbation theory, 
Casimir energy calculations provide a potentially powerful way to study nonperturbative 
effects. Several formally equivalent representations of the Casimir energy have been used 
in numerical calculations, notably (a) in terms of integrals over scattering phase shifts and 
bound states; (b) as an integral over the trace of a Greens function (equivalent to a functional 
determinant); and (c) as an infinite sum of Feynman diagrams. The aim of this paper is 
to study the various representations of the Casimir energy in a simple, highly convergent 
toy model, where the equivalence of different representations can be demonstrated without 
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the complication of divergences. Our analysis does not demonstrate the equivalence of 
representations in realistic models, where divergences make the arguments more complex. 
However, we do see how the different representations are related in a simple, calculable 
example. Also, because the model is so simple, it is possible to explore some issues of 
convergence and cancellation that are obscure in more realistic cases. 

The model we study makes use of the simple dynamics of separable potentials in nonrel- 
ativistic quantum mechanics. |I] Our first task is to explain how such a simple model can 
give insight into one-loop effects in a quantum field theory. Next we review the well-known 
Greens function and scattering theory analysis of a separable potential. With this in hand 
we discuss the various representations of the Casimir energy and explore their relationship. 
Finally, we study convergence of the perturbative expansion and associated computational 
issues. 

Formally, the Casimir energy is the sum of zero point energies for the modes of a quantum 
field, in the presence of some spatially varying background, (f>o(x), 

£ c = ±^(M0o]-^o) (1) 

where {u;[0o]} are the eigenfrequencies in the presence of the nontrivial background field 
0o, and {co> } are the eigenfrequencies in the trivial background, 4> = 0. The {co>[</> ]} are 
typically related to the eigenvalues of some simple differential operator, TC[(f>o], which looks 
like a Schrodinger operator for a scalar field, 

^="^+W (2) 

or a Dirac operator for a spinor field. V^, (x) is a "potential" derived from the field 4>q. 
Although we have written eq. (H) for a scalar in one- dimension, the same considerations 
apply in three-dimensions. Typically the {uj} are not themselves the eigenvalues, {A}, of 
7i but rather simple functions of the them. For example, in the case that \1/ is a scalar, 
ujj = a/ Xj + rri 2 . 

Our toy model is based on two alterations in this physical picture: first we assume that 
the {uj} are proportional to the eigenvalues of 7i, as they would be in the nonrelativistic 
case where ftw(k) = h 2 k 2 /2m (and k 2 is the eigenvalue of 7i); and second, we replace the 
local potential, V(x), by a nonlocal, but separable potential, 

V(x, x') = —Xu(x)u(x') (3) 

where J dxu 2 (x) = 1. Because of our first assumption the Casimir energy can be written as 



£ c = ±Tr(H-H ) 



(4) 
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As a consequence of the separability assumption, 7i can be written formally as 

H = H - X\u)(u\ (5) 

where \u) is the state with wavefunction u(x) = (x\u). Combining eqs. (||) and (0) we find 
that the Casimir energy is — |A, 

£ c = -±Tr{-A|u> (u\} 

= -|A . (6) 

This is the fundamental result that makes the study of the Casimir energy in separable 
potential models interesting - the answer is so transparent. Our object in this paper is to 
show that more conventional (and more complicated) methods of computing the Casimir 
energy coincide with this simple result. 

The remainder of the paper is organized as follows. In the next section we review 
the bound states, Greens functions and scattering amplitude for a separable potential. In 
Section [TTI| we derive expressions for the Casimir energy in terms of integrals over Greens 
functions or scattering phase shifts. In Section [IV] we compute the Casimir energy using the 
methods of Section |TTT| for the case of a separable potential and show that the result is — |A. 
In Section [V| we specialize to a particular choice of u and explore the convergence of the 
Casimir energy calculation as a function of the strength, A, of the interaction. We conclude 



in Section VI 



II. SEPARABLE POTENTIALS 

In this section we review the solution of the scattering problem for a scalar particle 
moving in an s-wave separable potential in three dimensions. We begin with the Schrodinger 
equation in three dimensions with a nonlocal potential, 

h 2 



2M 

where 



V 2 ^(r) + / d-W(r,r')^(r') = uip(r) (7) 



(r\V\r) = V(r,r) . (8) 



The resulting integral equation is in general more complicated than the local case, but for a 
separable potential, eq. (|3]), it is simpler. For simplicity we take a separable and spherically 
symmetric nonlocal potential, 

V(r, f') = -Xu{r)u{r') (9) 



R.L. JaSe and L.R. Williamson 



4 



for which eq. (0) reduces to 

h 2 



_ } . VV(r) - Xu(r) / dV u{r')ip{r') = uip(r) . (10) 

For A > the potential is attractive. 

A principal simplification with a spherically symmetric, separable potential is the absence 
of any interaction in partial waves with i > 0. This follows immediately from eq. ( |1~0"|) because 
J d 3 r' u(r')ip(r') projects out only the spherically symmetric part of ip. All partial waves 
except I = cancel out of the Casimir sum. We replace ip(r) by U(r)/r and find that 
eq. (|TUD simplifies to 



- U"(r) - An\ru(r) / dr'r'u{r')U(r') = k 2 U(r) (11) 
Jo 

where k 2 = 2Mu)/h 2 , and where we have set h = 2M = 1 henceforth. 

For scattering solutions, it is useful to convert eq. ( IID into an integral equation for the 
function, f/( + )(r), obeying scattering boundary conditions, 

POD POO 

C/(+)(r) = U (r) +4ttA / dr'G ( +) (r, r' , k)r'u{r') I dr"r"u{r")U {+) {r") (12) 



where Uo(r) is the free solution regular at the origin, Uq(t) = sin kr, and Gq ,(r, r', k) is the 
free, s-wave Greens function with outgoing wave boundary conditions, 

G { + \r,r',k) = ^e ikr > smk r< 

= Lr M j£*zL- (13) 

tt« J-oo (q ~k 2 - le) 

The properties of the Greens function guarantee that U^ + \r) satisfies the Schrodinger 
equation. The asymptotic form of the scattering wave at large r, 

lim U (+) (r) = -(e~ ikr - e 2im e lkr ) (14) 

! — >00 2 

and the behavior of the Greens function at large r enable us to read off the scattering 
amplitude in terms of U, 

f(k) = e i5(k) sinS(k) = —— \ drrsmkru(r) dr' r'U {+) (r')u(r') . (15) 



To proceed we return to eq. (|12"D, multiply by ru(r) and integrate over r. The resulting 
algebraic equation can be solved for J °° drru{r)U^ +S) {r), which yields f(k) upon substitution 
into eq. (|I~5D, 
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where we have defined 

/*oo poo 

£o(k) = / dr ru(r)U (r) = / drrsinkru(r) 
Jo Jo 

r°° 1 

X(co) = 8\ dq- -\Ul)\ 2 - (17) 

Jo q 2 -uj-te 

The scattering amplitude, f(k), has a cut along the positive k 2 axis induced by the cut 
in X(u). The cut begins with a branch point at threshold, k 2 = 0. The discontinuity across 
the cut is given by 

disc f(k) = 2i Im f(k) = 2i sin 2 5(k) (18) 

where we have used l/(x + is) = PV(l/x) — m5{x) to separate out the imaginary part 
when x is real. The phase shift S(k) can be read off eq. (0) conveniently by using the 
parameterization, f(k) = l/(cot5 — i), 

t «n 4w\Mk)\ 2 /k .... 

tan ^^ i-8A/ >| ?0 M|Vtf-^V (19) 

The Born approximation to the scattering amplitude is obtained by expanding the denomi- 
nator in eq. (|16D in a geometric series, 



n=0 

2 



q 2 — k 2 — ie 

In addition to the branch cut for real, positive k 2 , f(k) can be singular at values of k 



where the denominator in eq. (16) vanishes, i.e., where X(u) = 1. Bound states appear as 
poles in f{k) for k 2 < 0, or more precisely, k = in. [Poles at k — —in are "virtual states".] 
When k 2 < 0, X(— k 2 ) becomes real, 

\Uq)\ 2 

q 2 + k 1 



X(-, 2 ) = 8Xl dq^ 2 (21) 



and a bound state occurs at the value of for which 

= 8A r dg ^M! = l. (22) 
Jo y * o 

Since X(— k 2 ) is a decreasing function of k, the criterion for existence of a bound state is 
that X(0) > 1, or 

8A f°dq > 1 . (23) 

Jo 1 

This equation defines the critical value of A at which a bound state appears for a given 
choice of u(r). Note that there is at most one bound state in this separable potential. This 
completes our review of the properties of a particle moving in a separable s-wave potential. 
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III. THE CASIMIR ENERGY 

We are interested in (a) the Greens function and (b) the phase shift expressions for the 
Casimir energy. Both representations may be derived heuristically in the effective action 
formalism in field theory. [JJ] Equally well, we can start from the formal expression, eq. (|T]), 
and convert the sum over eigenenergies to a trace of a Greens function or an integral over 
phase shifts. 

A. Greens Function Representation 

The Greens function method starts from the Greens function in coordinate/energy rep- 
resentation, 

G{r,r ,u) = 2^— — — ( 24 ) 

^ UJj — UJ 

j J 

where the {(fiji?)} are the unit-normalizedQ energy eigenstates in coordinate space and the 
summation ranges over the spectrum of TC, so 

(H - uj)G{t, r', lu) = 5 3 {r - f') . (25) 

In order to relate this to the Casimir energy, consider the difference of the Greens functions 
evaluated at uo — it and uj + ie integrated over space, 

d 3 r [AG(r, r,u + ie) - AG(r, r, uj - ie)} = 2ixi Y^[5(^ - uj) - 5(uj - co 0j )] (26) 

3 

where AG is the difference between the Greens function in the background potential and 
the free Greens function. To obtain the Casimir energy, multiply eq. (B^) by uj, integrate 
from — oo to oo, and divide by 4m, 

-1 /'OO /"OO 

£ c = — dujuj d 3 r[AG(r,r',uj + ie)- AG(r,r',uj-ie)] . (27) 
J_ 00 Jq 

This result can be simplified further by introducing a function, F(uj), defined by its 
derivative, 



dF(uj) 
diu 



d 3 rAG(r,r,uj) (28) 



We proceed as if the spectrum is discrete and the states are normalizable. Our result applies 
equally to the case (of interest) where the spectrum is continuous and the states are normalized to 
5- functions. 
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FIG. 1. Contours in the complex w-plane relevant to evaluating the Casimir energy. 



and the condition that F(u) 
ing by parts, we find 



as \u\ — ► oo. Then substituting into eq. (|27|) and integrat 



Sr. 



1 

Am 



duj [F(u 



F(u + ie)} 



(29) 



The surface terms at ±oo ± ie can be shown to cancel. 

Referring back to the definition of AG, we see that the discontinuity in F along the Re 
to axis is associated with the eigenstates of 7i. Let u be the ground state of TC. For all 
u) < u we can send e — > 0. Then the integral of eq. ( p9|) can be written as a contour integral 
over a contour C (as shown in Fig. |J) that runs in from +oo above the Re u axis to u>o and 
then returns to +oo below the Re uo axis, 



6 Ani 



duF(u) — ——7 (t duoFiuj) . 
Am 1 



(30) 



In the second term we have replaced the contour C by the counterclockwise circle at infin- 
ity, Cqo, since in general the only singularities in F(u) lie on the real axis. This compact 
expression will be particularly useful in the separable case. 



B. Phase Shift Representation 

The phase shift representation for £ c begins from eq. (||) divided up into bound state 
and continuum contributions, 



dnuj 

' continuum 



bound Jconth 
states 

poo 

= ~\y2 B 3 + \ I dup(uj)uj 
Jo 



(31) 
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where Bj are the binding energies, and p(uj) = dn/duj is the modification of the continuum 
density of states due to the interaction. A simple and general argument connects p{yj) to 
the phase shift (in our case only the s-wave contributes), 

. , 1 d5(uj) . . 

P w = --^ . 32 
7r do; 

Substituting p(u;), and integrating by parts yields, 

1 f°° 

£ c = -\H B l-^ J o dUjS ^ ■ ( 33 ) 

The surface term in the integration by parts ~ uj5(uj) vanishes trivially at uj = 0. For a 
general problem the surface term at infinity causes problems because generically o~{uj) ~ 1/uj 
as uj — > oo. However, in the case of a separable potential, we shall see that the phase shift 
vanishes rapidly as uj — > oo, so the surface term can be ignored. So we may immediately 
calculate the Casimir energy if we are given the phase shift. Eq. (|3~0D and eq. ([33|) give 
alternative representations of the Casimir energy, which should both be equivalent to eq. © 
for the separable case. 



IV. CASIMIR ENERGY FOR THE SEPARABLE CASE 



In this section we examine the Green's function and phase shift representations of the 
Casimir energy for the separable case. 



A. Greens Function Representation 

To exploit eq. ( |2"7| ) we need an explicit expression for the trace of the Greens function 
for the separable potential of Section 2. We start with the Lippmann-Schwinger equation 
for G, 

G(r,r',uj) = G (r,r',uj) + A J d 3 r t d 3 r2Go(r,r 1 ,uj)u(r 1 )u(r 2 )G(r2,r ,uj) (34) 

which may be solved by iteration, 

AG(r,r',uj) = \ J d 3 n d 3 r 2 G (r, r u uj)u{r) [1+X{uj)+X{uj) 2 + ■ ■ ■] u(r')G (r 2 , r', uj) (35) 

where X(uj) was defined in eq. (|17|). 

Eq. (|27D requires the us to set f = r' and integrate. If we substitute from eq. (|35|) the 
result takes a particularly simple form, 
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J d 3 r AG(f,r,u) = ^ ^2 



dx ^ [x(u)Y 



m 

m=0 



-A-ln(l-X(u)) (36) 



so 



S c = <b duo\n(l-X(uo)) (37) 

where we have identified the generic function F(u) of eq. ( p8|) with — In (1 — X(u)) in the 
separable case. [With this identification it is easy to check that the surface terms discussed 
following eq. (|2"9]) do indeed vanish.] To evaluate the integral we need lim^oo In (1 — X(u)). 
From eq. ( |T7D and the normalization of u(r), it is easy to see that asw^ oo, 

In (1 - X(u)) ~ -X(uo) ~ A/cj . (38) 

Thus, 

S c = & — = — . 39 

Am J Coo oo 2 K ' 

Confirming the value of the Casimir energy is —A/2 as we derived in Section 1 from more 
formal considerations. 

Some comments are in order: 

• The method we have presented is equivalent to the usual, graphical analysis of the 
effective action in quantum field theory. Had we begun with a theory describing a 
scalar field, propagating in a scalar background, $, then the Casimir energy would 
have been given by the (infinite) sum of one-loop graphs for \1/ with insertions of 
the coupling. The one- loop graphs could then be resummed into a Lippmann 
Schwinger equation for G(r, r', uo). The special advantages of a separable potential are 
(a) that the resulting integral equation for G is solvable, and (b) that there are no 
ultraviolet divergences (seen as divergences in the uo — ► oo limit) to complicate the 
calculation. 

• The resummation of 1 + X + X 2 + ■ • • that generated the logarithm in eq. ( |3"5| ) is valid 
for small A where the series converges. The result can then be analytically continued 
to large A were the re-expansion into a (Born) series does not converge. 



B. Phase Shift Representation 



In Section 3 we also derived a simple representation, eq. fl3"3]), for the Casimir energy as 
a sum over binding energies and an integral over the phase shift 8 (id). In this subsection we 
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show explicitly that this representation is equivalent to the Green's function representation 
as found, for example, in eq. fl3~7|). 

First, we return to the contour, C, of Fig. [I], 

£ c = — — Idu ln(l -X(u)). (40) 
4tu Jc 

The contour integral is equivalent to integrating the discontinuity in the integrand across 
the real axis from just below the lowest bound state uj q to oo, 

i r°° 

g c = — / da;Imm(l -X(u + ie)) (41) 

where we have used the fact that the discontinuity in the logarithm is 2i times its imaginary 
part. 

There are two distinct regions in the integral. For u < 0, X(u) is explicitly real, and 
the logarithm can have an imaginary part if and only if X(u) > 1. For u > 0, X(u) always 
has an imaginary part (see eq. (|T7|)) related to the scattering amplitude. We shall show 
that these two contributions map into the binding energy and integral over the phase shift 
respectively, as expected from eq. fl33|) . 

First consider uj < 0. According to the analysis of Section 2, eqs. (|2~l[) and (|22), etc., 
X(u) takes its maximum at uj = 0, and decreases as ui decreases. If -X"(0) > 1 there is a 
single bound state with binding energy, B = —ujq, determined by the equation X(—B) = 1. 
Thus if A is large enough to generate a bound state, then In (1 — X(u)) has an imaginary 
part (equal to it) for — B < uj < 0. Thus 

1 r° 

— (kulm\n(l-X(uj)) = -\B . (42) 
Next consider u > 0. From the analysis of Section 2, we find 

In (1 - XM) = In (l - 4A-T dq JM^- + 4^^^ . (43) 
From this we read off the imaginary part, 

which is just — tan8(k) as defined in eq. (^9[) . Substituting from eqs. ( fl2"D and (|44] ) into 
eq. (|4TD , we confirm the phase shift plus binding energy representation, eq. ([53]) , 



1 f 00 

^ = / ^^)- ( 45 ) 



2tt 

3 
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V. CONVERGENCE AND NUMERICAL ISSUES IN A SEPARABLE 

POTENTIAL MODEL 

A particular virtue of the separable potential model is that it is simple enough to allow 
us to investigate issues that are difficult to attack in more realistic theories. As an example 
of such an issue, we consider here some computational aspects of phase shift representation 
of the Casimir energy. Ref. introduced a method of computing the Casimir energy 
in realistic quantum field theories based on the Born approximation to the phase shift. 
Divergences generated by the lowest-order Feynman diagrams are associated with the first 
few terms in the Born approximation to the phase shift. To remove the divergences from 
the numerical part of the calculation, the first few Born approximations are subtracted from 
the phase shift. These contributions are then added back into the calculation in the form of 
(divergent) Feynman diagrams, which are regulated and renormalized by means of available 
counterterms. Schematically 



where 8 is the phase shift with the first N terms in the Born approximation subtracted, 



D n is the contribution to S c from the n -order Feynman diagram, and CT denotes the 
renormalization counterterms. Both the {D n } and the counterterms depend on a cutoff if 
the theory has ultraviolet divergences. The Born-subtracted phase shift, however, is cutoff 
independent. 

In realistic quantum field theory applications the Feynman diagrams, {D n }, and the 
counterterms are calculated by standard, diagrammatic methods. The integral over 8 N 
is done numerically. The Born approximations {X n 8^} have just the right behavior to 
cancel the large uo tail of 8(u) and render the tu-integral convergent. However they are a 
poor approximation to 8 (us) for small u, especially in partial waves where there are bound 
states. P| Numerical studies show that at low u the function 8 (u) is much larger than 
its integral J °° dw 8 N {ui). The origin of this effect is obscured in those studies because 
the integral over the individual Born approximations, J °° du <r n -*(u;), diverge. The large 
uj behavior of the separable potential model is very soft, enabling us to study this issue 
explicitly. 

It is convenient to specialize further to a particular choice of the separable potential 
function, u(r). Following Gottfried, we choose the Yukawa function, 




(46) 



N 




(47) 



n=l 
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„(r) = . (48) 



A simple calculation gives 




a k 



2tt k 2 + a 



2 



• Y « = 2{^w ■ (49) 

It is convenient to use the momentum, k, rather than u = k 2 as the independent variable, 
also to introduce scaled variables, k/a = z, X/a 2 = g, and finally, it is easiest to explore the 
model by studying the scattering amplitude, f(k) — > f(g, z) = 1/ (cot S(g, z) — i), 

(1 + iz) z ((l — iz) 1 — g) 

The scattering amplitude f(g,z) has two types of singularities. The double pole at 
z = —i is a "potential singularity" arising from the fourier transform of u(r), £o(&)- It is 
on the second sheet of the complex u plane where it does not affect the contour shifting 
arguments we used earlier in this paper. The other singularities are poles at 



z ± ={ • (51) 




For g > 1, z+ is associated with a bound state: it lies on the positive imaginary axis, 
corresponding to the negative real axis on the first sheet of the complex k>plane.f] The 
energy of the bound state is ooo = —a 2 z 2 _ = — (y/X — a) 2 . The pole at z_ is always on the 
second sheet and is not dynamically important. The singularities in the complex-z plane 
are summarized (for g > 1) in Fig. ^. 

With these considerations in mind, we can write an illuminating expression for the 
Casimir energy, 

{--!- J duo5(uo) for A < a 2 

^ 1 poo (52) 

-\(V\-af - — J duj5(uj) for A > a 2 

i.e., where there is a bound state, it appears explicitly in S c . The Born approximation is an 
expansion of S(g,z) in powers of g obtained by expanding the denominator of f(g,z) (see 
eq. (|5C)D ) in a geometric series: 



2 For g < 1 the z+ pole lies on the negative imaginary axis, i.e., on the second sheet of the w-plane. 
As g — > 1 it becomes virtual state on its way to becoming a bound state. 
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FIG. 2. The Born expansion converges outside a circle of radius yfg centered at —i in the complex 
z-plane. Singularities in f(g,z) at z = —i, z + = i(^/g — 1), and z_ = — i(y/g + 1) are shown. 



(53) 



The convergence of the Born expansion is determined by the convergence of the geometric 
series for f(g,z), which converges provided |1 — iz\ > y/g. That is, z must lie outside the 
circle of radius y/g centered at z = —i. To compute the Casimir energy by the phase shift 
method, we must integrate over all u > 0, corresponding to z G [0, oo]. 

As long as g < 1 (A < a 2 ) the Born expansion converges for all real, positive z and can 
be used under the integral sign in eq. fl52|) to compute S c , 



1 00 PC 



dtu5 {n h 



LU 



for A < a z 



(54) 



This is a remarkable result: both sides are power series in A, but the left hand side consists 
of a single term. We conclude that, for A < a 2 , when there is no bound state, the entire 
Casimir energy is given by the first Born approximation, and that the integrals over all 
higher Born approximations to the phase shift vanish: 

A 



Z7T 



UJ) 



dtoS (n \ 



UJ 



for all n > 1 



(55) 
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FIG. 3. The phase shift for the Yukawa separable potential for g = 0.99. The solid curve is the 
exact result, the four dashed curves give the first four Born approximations respectively. 



This result may easily be verified numerically. 

When g > 1 (A > a 2 ), it is no longer possible to use the Born approximation to evaluate 
the cu-integral. A glance at the second equation in eq. (|52|) confirms that the Born approxi- 
mation cannot converge when there is a bound state: if it did, it would yield a result analytic 
in A, however the cu-integral must generate a factor a^/\ to cancel the contribution of the 
bound state and give S c = — |A. To generate V\, the Born expansion must not converge. 
Still, it is useful to subtract some number of terms in the Born expansion from the exact 



phase shift (5 



N i 



S(co) — YlnLi X n 5^ nS) (to)) and to rewrite £ c accordingly, 



Sr. 



^ Jo 27r Jo 



duj 5 



N, 



1 

2^ 



dou 5 



N / 



(56) 



Here we have subtracted the first iV terms in the Born approximation to convert 5 to 
5 N . When we added them back in, we used the fact that the integral of the first Born 
approximation gives exactly — |A, and that the integrals of all higher Born approximations 
vanish. From this analysis we conclude that when there is a bound state, the integral of the 
"Born subtracted" phase shift, 5 N gives an iV-independent contribution to S c , 



1 (L:!) 



a) 



(57) 



These results help to understand the effects that have been seen in numerical calculations 
in more complicated theories. In our toy model, when g is such that there is no bound state, 
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FIG. 4. The phase shift for the Yukawa separable potential for g = 4.5. The solid curve is the 
exact result, the four dashed curves give the first four Born approximations respectively. 



the Casimir energy is given entirely by the first Born approximation (which, in a theory 
with divergences, would be computed from a Feynman graph). The n th Born approximation 
gives a contribution to 5 of order g n , which may be large for g < 1, but its integrated 
contribution to the Casimir energy is exactly zero. In more complicated theories the answer 
is very small but not exactly zero. The magnitude of the effect can be seen in Fig. |3], 
where we plot the exact phase shift and the first four Born approximations for g = 0.99. 
The first Born approximation has the same integral as the full phase shift. Each higher Born 
approximation integrates to zero. The situation becomes even more bizarre in when there is 
a bound state. The entire Casimir energy is still given by the first Born approximation! The 
bound state contribution is exactly canceled by the remainder of the phase shift integral. 
One can remove any number of further Born approximations from 8 - defining 8 , which is 
significantly modified by the subtraction of higher Born approximations - without changing 
the integral. The effect can be seen in Fig. where we plot the exact phase shift and the 
first four Born approximations for g = 4.5. Once again the toy model caricatures effects 
seen in more complex and more realistic theories. 



VI. SUMMARY AND CONCLUSIONS 



Our principal purpose has been to develop a simple model to obtain insight into com- 
plicated calculations of Casimir energies in quantum field theories. We have verified in our 
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separable potential model that several traditional representations of the Casimir energy em- 
ployed in quantum field theories are equivalent. Our model is very simple. In particular, it 
does not display the divergences characteristic of real quantum field theories. Thus we are 
not surprised that results are confirmed in our model that might be spoiled by divergences 
in more realistic theories. Nevertheless, the model is rich enough to illustrate some of the 
computational difficulties that have been seen numerically in more complicated theories. 

Also, we should point out that there are cases where separable potentials do arise in the 
treatment of more realistic theories. In particular, Bashinsky has shown that the collective 
coordinate quantization in theories with solitons introduces additional, separable-potential- 
like terms into the small oscillations Hamiltonian whose eigenvalues determine the Casimir 
energy. Q The ideas developed here would apply relatively directly to such cases. 
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